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X) ■ ABSTRACT 

Context. Third-order galaxy-galaxy lensing (G3L) is a next generation galaxy-galaxy lensing technique that either measures the 

■ excess shear about lens pairs or the excess shear-shear correlations about lenses. From their definition it is clear that these statistics 
' assess the three-point correlations between galaxy positions and projected matter density. 

Aims. For future applications of these novel statistics, we aim at a more intuitive understanding of G3L to isolate the main features 
that possibly can be measured. 

Methods. We construct a toy model ("isolated lens model"; ILM) for the distribution of galaxies and associated matter to determine 
the measured quantities of the two G3L correlation functions and traditional galaxy-galaxy lensing (GGL) in a simplified context. 
(-H \ The ILM presumes single lens galaxies to be embedded inside arbitrary matter haloes that, however, are statistically independent 

^ J ("isolated") from any other halo or lens position. Clusters of galaxies and their common cluster matter haloes are a consequence of 

I ■ clustering smaller haloes. In particular, the average mass-to-galaxy number ratio of clusters of any size cannot change in the ILM. 

O ' Results. GGL and galaxy clustering alone cannot distinguish an ILM from any more complex scenario. The lens-lens-shear correlator 

^ I in combination with second-order statistics enables us to detect deviations from a ILM, though. This can be quantified by a difference 

C/5 . signal defined in the paper. We demonstrate with the ILM that this correlator picks up the excess matter distribution about galaxy pairs 

■ inside clusters, whereas pairs with lenses well separated in redshift only suppress the overall amplitude of the correlator. The amplitude 
suppression can be calibrated. The shear-shear-lens correlator is sensitive to variations among matter haloes. In principle, it could be 
devised to constrain the ellipticities of haloes, without the need for luminous tracers, or maybe even random halo substructure. 

^ ' Key words. Gravitational lensing:weak - Galaxies: halos - (Cosmology:) large-scale structure of Universe 

(N 

'^1. Introduction 

■ Gravitational lensing (Schneider 2006, for a recent review) has established itself as valuable tool for cosmology to investigate the 

■ large-scale distribution of matter and its relation to visible tracers such as galaxies. In the currently favoured standard model of 
Ps) , cosmology (e.g., Peebles 1993; Dodelson 2003), the major fraction of matter is a non-baryonic cold dark matter component, i.e., 
T—i • matter with non-relativistic velocities during the epoch of cosmic structure formation. In this situation, gravitational lensing is an 

^ I excellent probe as it is sensitive to all matter as long as it interacts gravitationally. 
• I— I ■ The main observable in lensing is the distortion of shapes of galaxy images, in the weak lensing regime mainly "shear", by the 

■ intervening inhomogeneous gravitational potential that is traversed by light bundles from the galaxy. Over the course of the past 
J-j ] decade, applications of the gravitational lensing effect have come of age. To name a few results (Bartelmann 2010, and references 

therein), it was used to map the dark matter distribution, to study the matter density profiles in galaxy clusters and to determine their 
masses, to measure the relation between the galaxy and dark matter distribution (the so-called galaxy bias), to constrain the total 
matter density of the Universe and its fluctuation power spectrum, and very recently to gather independent evidence for the overall 
accelerated expansion of the cosmos (Schrabback et al. 2010). 

Of particular interest for this paper is the so-called galaxy-galaxy lensing technique where positions of foreground galaxies 
("lenses") are correlated with the (weak lensing) shear on background galaxy images ("sources"). Thereby, statistical information 
on the projected matter distribution around lens galaxies can be extracted. Since the first attempt by Tyson et al. (1984) and the 
first detection (Brainerd et al. 1996; Griffiths et al. 1996) of this effect, galaxy-galaxy lensing nowadays is a widely applied robust 
method to study the galaxy-matter connection (Fischer et al. 2000; McKay et al. 2001; Guzik & Seljak 2002; Hoekstra et al. 2002; 
Pen et al. 2003; Hoekstra et al. 2004; Seljak & Warren 2004; Sheldon et al. 2004; Mandelbaum et al. 2005; Kleinheinrich et al. 
2005; Mandelbaum et al. 2006b, a; Simon et al. 2007; Parker et al. 2007; van Uitert et al. 201 1). The traditional and hitherto mainly 
employed approach is to correlate the position of one lens with the shear of one source galaxy (GGL hereafter). 

Schneider & Watts (2005) advanced the traditional GGL by considering a new set of three-point correlation functions that 
either involves two lenses and one source ("correlator 0") or two sources and one lens ("correlator G±"). Both correlators are tools 
to directly study higher-order correlations between galaxies and the surrounding matter field. In the literature, this technique is 
termed 3"'-order galaxy-galaxy lensing or galaxy-galaxy-galaxy lensing. There are alternative but mathematically equivalent ways 
to express these statistics, e.g., the aperture statistics (N^M^p) and {NM^^} instead of the three-point correlation functions G± and 
(Schneider & Watts 2005; Simon et al. 2008). In practical measurements usually the aperture statistics are preferred, as they 
automatically remove unconnected 2"''-order contributions in estimators of the statistics and allow one to separate E-modes from 
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B-modes or parity modes, of which the latter two cannot be generated by gravitational lensing as leading order effect. In this paper, 
we focus on the E-modes in the correlators Q and G±. 

G3L has already been measured in contemporary lensing surveys, such as the Red-Sequence Cluster Survey (Simon et al. 2008), 
and will therefore presumably be routinely measured with ongoing surveys such as KiDS', Pan-STARRS^, DES^, or in the future 
surveys Euclid"* and LSST''. The prospects of learning more on the galaxy-matter relation or new observational tests for theoretical 
galaxy models (e.g. Weinberg et al. 2004; Bower et al. 2006; De Lucia et al. 2007) with G3L are thus quite promising. 

The main obstacle for exploiting the new G3L statistics is their physical interpretation. It is clear from the definition that Q 
quantifies the shear signal (or projected matter density) in excess of purely randomly distributed lenses, picking up only signal from 
clustered lens pairs (e.g., Johnston 2006), and that G± is a two-point correlation function of shear associated with matter physically 
close to lenses. It is unclear, however, what physical information this translates to and what new feature may be contained in G3L 
that may be missing or is degenerate in traditional GGL. To elucidate these new statistics and to pave the way for new applications of 
G3L, we conceive here a simplistic model for the distribution of lenses and matter: the isolated lens model. Then G3L is flashed out 
in the light of this model. For the definition of quantities relevant for weak gravitational lensing, we refer the reader to Bartelmann 
& Schneider (2001). The mathematical machinery of a halo model expansion devised in the calculations is very similar to Scherrer 
Sl Bertschinger (1991), although used in a different physical context. 

The structure of the paper lays out as follows. Sect. 2 introduces our model and derives the tangential shear about a lens, the GGL 
signal, expected from this description. Sect. 3 moves on to calculate the lens-lens-shear or Q correlator for this specific scenario. 
Sect. 4 does the same for the shear-shear-lens or G± correlator The final Sect. 5 summarises the main conclusions drawn in the 
preceding sections. 



2. The isolated lens model and galaxy-galaxy lensing 

Here we lay out a simple model for the distribution of matter and galaxies inside it. This model is founded on the assumption that 
lenses are embedded inside a matter halo that generates the shear profile yb{0\ a) acting upon a background source, where is the 
separation vector from the centroid of the galaxy, and a denotes a set of intrinsic halo parameters that control the matter density 
profile of the halo. Importantly, the intrinsic parameters are statistically independent of the intrinsic halo parameters of any other 
halo or the separations of other lenses. We hence coin lenses and their host haloes in this scenario "isolated". In this sense, this is a 
very crude halo model representation (Cooray & Sheth 2002) of the lens and matter distribution, assuming for simplicity that every 
halo is occupied by exactly one (lens) galaxy. Notice that matter which is statistically independent of the lenses does not need to be 
accounted for, as this would not contribute to a galaxy-matter cross-correlation function, although it certainly would affect the noise 
in a measurement. 

In the following, 2D positions on the flat sky are, for convenience, denoted by complex numbers - 6i + iOj where 0\ {62) is 
the position in direction of the A:()')-axis. By = \0\ - y/00* we denote the modulus of 0. Likewise the Cartesian shear 2-spinor 
Jc - 7i + iy2 is denoted as complex number. We define the tangential, ji, and cross, yx- shear of yc(^) relative to the origin by 

y(0; if) := y,(0) + ij^(0) = -^-^"^7^ , (1) 
where <f is the polar angle of 0. 

2.1. Isolated lens model 

Within the isolated lens model (ILM hereafter), the number density distribution of lenses on the (flat) sky with area A is 

«gW = Xi^D(^"^5')> (2) 

1=1 

where 6^^\0) is the Dirac delta function, and the resulting shear field is 

rc(e) = 2rh(0-^;«,), (3) 

!=1 

Sticking a shear profile to every of the A^d lens position (f}. The shear profile is directly related to the projected matter density about 
the lens. In the following we will use the lens number density contrast 

K.{0) - 1 , (4) 



' http : //www . astro-wise . org/pro jects/KIDS/ 
^ http://www.cfht. Hawaii .edu/Science/CFHLS/ 
' http://www.darkenergysurvey.org 

http://sci.esa.int/euclid, see also Laureijs et al. (2011) 
' http://www.lsst.org/ 
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Fig. 1. Illustration of the isolated lens model. Galaxies are depicted by black pixels, their matter haloes as red disks. For simplicity all haloes 
are identical in this visualisation. Left panel: Lenses are distributed randomly on the sky. Right panel: Lenses cluster to produce clumps with a 
common matter envelope, yet still described as sums of individual haloes. On the statistical level, this clumping is quantified by non-vanishing 
2"''-order and 3'''-order correlation functions oj and £1, respectively. 

where «g := Na/A is the mean number density of lenses within the area A. For the 2"''-order angular clustering correlation function 
of lenses on the sky (e.g. Peebles 1980), we employ the function 

w(l^l2l) = (^g(eiM^2)) (5) 

with 0ij :- Oi - 0j being the separation vector of two positions. A value uiOn) > expresses an excess of galaxy pairs at separation 
012 compared to a purely random distribution. 

The ILM is a more general description than may appear at first sight: Galaxy clumps with joint matter envelopes ("galaxy 
clusters") are not explicitly excluded, although every individual halo does host only one lens. To form matter haloes of whole 
clusters, we can always stick together and overlap matter haloes, changing the clustering correlation functions in consequence. The 
difference between the left and right panel of Fig. 1 lies therefore in the choice of the clustering correlation functions, which will 
enter the following calculations. Hence, the ILM expands larger matter haloes as sums of individual matter haloes of clustered 
galaxies. Crucially, however, the ILM is incapable to implement clumps of galaxies that contain on average more mass than 
isolated galaxies; the mean matter-to-lens-number ratio has to be constant throughout. To form a clump with a higher mass-to-lens 
number ratio would require to increase the mass of all individual haloes simultaneously, i.e., to change the internal halo parameters 
of all lenses inside the clump in a similar fashion. This is not allowed, except by chance, however, since internal parameters are 
statistically independent. A generalisation of the ILM in this direction could be achieved in a full-scale halo model with matter 
haloes hosting more than one galaxy. Note that the luminosity of a lens can also be seen as an internal parameter Therefore, a 
constant mass-to-lens number ratio plus statistical independence of internal parameters amounts also to a constant mass-to-light 
ratio of all clumps in the model. 



2.2. Galaxy-galaxy lensing 

Before we embark on 3"^-order statistics, we start with the more familiar 2"'^-order GGL (Bartelmann & Schneider 2001). These 
statistics are the mean tangential shear about a lens at separation § - \02 - 6\\ (Fig. 2), defined by the correlator 



«g(0l)re(02) 



.+2i^j 



=0 



(6) 



Owing to isotropy, the mean tangential shear is only a function of separation i? and independent of the polar angle (f. The bracket 
{. . .) denotes the ensemble average over lens number densities and shear configurations. The underbraced term has to vanish due to 
the statistical isotropy and homogeneity. 

As lens number densities and shear configurations are expanded in terms of haloes in the ILM, we consider the correlator as 
ensemble average over aU possible lens positions and internal halo parameters: 



e))y^{e2-e)-aj) 



■ ^)y^{e: 



2 - 0); ff/)) + Yj (-^d'^^i - ^')rh(»2 - e)- aj)j 



(7) 
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Fig. 2. Left panel: Illustration of the parametrisation of the GGL correlator 7,(1?). The lens is located at 61, the source is at 62- Middle and 
right panel: Illustration of the parametrisation of the G3L three-point correlators G+(&i , ^3) (middle), and the galaxy-galaxy-shear correlation, 
0i&i,&2, (pi) (right). The figures are copied from Schneider & Watts (2005). 



which splits into two separate sums with ensemble averages over all halo parameters of one halo (one-halo term) 
or two haloes (two-halo term) 

The statistical independence of halo positions 0^ and internal halo parameters or, is explicitly used here; (. . and (. . ■)ai,aj the 
averages over internal halo parameters of a single halo or jointly for two haloes, respectively. For the latter, we stress again that we 
will assume statistical independence of or, and aj. 

In the following, we will need the average halo shear profile 



(8) 



(9) 



rhW:=(rh(e;«))^ = j daP„(a)rh(0;a) 



(10) 



Paia) is the probabihty density distribution function (p.d.f.) of the internal halo parameters a. The tangential halo shear profile y,jj 
is not to be confused with y^ in Eq. (6) that describes the total mean tangential shear about a lens including contributions from the 
lens halo and haloes of clustering neighbouring lenses. Due to rotational symmetry, the average profile 7h h^s a vanishing cross 
shear component, for which reason we can express it in terms of the tangential halo shear function ^{6), which is only a function 
of the separation Q. 

Utilising this definition, we arrive for the GGL correlator at; 



(11) 



=0 



The underbraced term must vanish due to radial symmetry. Inside the sums all terms become independent of individual lens positions 
9^ and halo parameters a, due to the averaging. The last step assumes that the number of haloes is large, i.e., A^d » 1, in particular 
Nd{Nd - 1) ~ N^- This will also be assumed for all following calculations. 
Employing (6) and (10), we finally find: 



f r>oo /-'2n 

= rth('?) + «g d0'0'd^'e-*'^'-^'w(>I')r,,((?') = rth('?) + «g d06y,^(0) d^cos(2.^)6>(>10 
J Jo Jo 



y,(§) 
with the expression 



(12) 



(13) 



The last step in (12) exploits that cl){6) has vanishing imaginary part. In the specific ILM description, the average shear is expanded 
in terms of two components: The first term in Eq. (12) is the one-halo term, dominating at small separations, whereas the second 
term is the two-halo term due to the clustering of haloes. Importantly, GGL is only sensitive to the average lens halo yj,(^) but 
insensitive to deviations of yh{0; a) from y^{9) in the actual halo population, which are explicitly allowed within the ILM. 
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3. Lens-lens-shear correlator 

We now turn to the 3"'-order galaxy-galaxy lensing statistics, starting with the constellation of two lenses and one source as depicted 
in the right panel of Fig. 2. The correlator considers a cross-correlation between lens number densities at two positions 0i and 02 
and the shear at 6^ : 

lnJ0i)nJ02)ydO3)) = -«;e+'<'^'+^=' ^2, ^3) ■ (14) 



Statistical homogeneity and isotropy implies that we can extract a correlation function ^ from the correlator that is solely a function 
of lens and source separations, §1 and §2, and the opening angle ip^. This function is split into an unconnected, Qnc, and a connected 
part 

gi&U&2,h) = (^g(0lK(02)r(^3;^^^)) =^(>?l,'?2,03)-^nc('?l,'?2,^3) - (15) 

The connected part vanishes for unclustered lenses, while the unconnected part can be shown to be generally (Schneider & Watts 
2005): 

&nMx,-&2,4>3) ^-'"^'%m + &^'"^'%m , (16) 

which is just the sum of the GGL shear profile around each lens. Therefore, Q encodes the shear in excess of what is expected from 
unclustered lenses with the average shear profile around them, or: It quantifies the shear signal about clustered lens pairs. Note 
that the unconnected terms do not contribute to the aperture statistics (N^M^p) and are thus not the primary quantity measured with 
G3L. 



3.1. Derivation 

For within the ILM, we need to evaluate the connected terms of the correlator («g(0i)«g(02)Tc(^3)) with the model specifics Eqs. 
(2) and (3): 



(n^(0i)n^(02)7c(03)) = ( £ ^D^(^i - «?)^d'(^2 - ^)rh(«3 - ^; «a)) (17) 
= (^d*(^i - «?)'^d'(^2 - e^)jm - ^\ ad) + te'(«i - »')^d'(«2 - e))yu(03 - ; ff,)). . + C2 perm. ) 



one-halo two-halo 



where 

{■■■),.,■■- J3 / d'(^dW(^k [ 1 + ^((^j) + ^(Sl) + ^(g).) (^^)][ ...)^ ^ ^ (18) 

unconnected 

is the ensemble average over three haloes (three-halo term). By 

^(013, 023, 012) = (^g(0l)'Cg(02)^g(03)) (19) 

we denote the (connected) 3'^''-order angular clustering correlation function of the lenses that only depends on relative galaxy 
separations. The sum (17) hence decays into a one-halo term (first sum), two-halo (next sum plus two identical sums apart from 
permutations of the indices / and j) and the three-halo term (last sum). The one-halo term vanishes for 0i + 02, i.e., distinct lens 
positions, owing to the Delta functions. For the connected terms, in the halo correlators (. . or (. . .),j;jt only the summands with 
the leading order clustering correlation functions are relevant, i.e., only the terms with Q in the correlator (. . .),,/,*:, while terms 
with CO or 1 are part of the unconnected terms (underbraced). They belong to the unconnected part of the 3'^''-order lens clustering. 
Similarly, in {. . .), y only terms associated with o) are of relevance, the terms generated by 1, the unconnected part of the 2"''-order 
lens clustering, will go into Q^^^. 

By evaluation of all these ensemble averages one thereby obtains for 0i + 02 + 03'. 

G{M2,<p3)^ §2^{-»l , )?2, 03) + ^3h(l?l , 1?2, 03) (20) 
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with the two-halo terms (we utilise the relation Jhi-O) - 7'h(^) following from Eq. 10) 

- a;(i?3)(e-'^7,,h('?i) + e^'^'rt.h('?2)) , (21) 
the lens-lens separation 



&3 = ^&-^+d-l-2§i§2 COS 03 (22) 

and the three-halo term 

03h(^u^2, 03) := -nge-'^'^'^'^^' j (f0D.{\0i - 0\, \02 - 0\, \0i - 02\)%{03 - 0) 
= -Hge-'<'^'+'^^' J d^6n{\0n+0\,\023 + O\,&3)n(O) 

= n„J d6'6id,^Q(Y(i?i,0,<^-^Oi),Y(i?2,6i,</'-^2),)?3)e"'<^'+^^"'^>r,,i,(f) 
= «g J d00d^Q(Y(i?i,0,^ + 03),m,0,¥'X'^3)e+-'*'+^>rt,h(0) 

d6'6ir,h(6') j d^cos(2^)Q(T()?i,6»,^ + 03),T()?2,6i,^),j?3) , (23) 
for which we have introduced the auxiliary function 

Y(6'i,6'2,0) := 6^+ el +20102 cos (f> ■ (24) 

The transformations in (23) use <pj, - ip2-if\ and a change of the integral variables ^ i-> (f + (f2 and i-> + 0^. The last step utilises 
that . ■) has a vanishing imaginary part. 

3.2. Interpretation 

In the context of 0, we can define a exceii /izoii map in the following way. The function Q can, for lens pairs of fixed separation i?3, 
be mapped as excess shear field at position ^3 with Cartesian shear value 

70(^31^1,^2) - , ^(013,023.03) , (25) 

|t'13l|i'23l 

where <pT, is defined as angle spanned by 023 and ^13. In this map, we fix the lens positions at 0i = +i?3/2 and O2 - -t?3/2 on the 
X-axis. This shear map can be converted into a convergence map (e.g. Kaiser & Squires 1993), as for example done in Simon et al. 
(2008). The excess mass from the two-halo term of the ILM is just the halo mass about each lens at 0i and 02, weighed with the 
clustering strength (y(i?3) of the pair This is exactly the mass we would anticipate around a pair of lenses, after one has subtracted 
the mass around unclustered pairs (Eq. 21 with <y(i?3) si) and if one ignored the effect of third haloes. The excess mass stemming 
from a third halo, clustering around the lens pair, is described by the three-halo term. 

That the excess shear or mass originates from galaxy clusters can be argued from the ILM (Fig. 3). In the ILM, the excess 
shear is expressed in terms of the 2"''- and 3"'-order lens clustering correlation functions on the sky. We have no signal, if lenses are 
unclustered, i.e., w = Q = 0, or simply if we have no galaxy clusters. Unclustered lenses with statistically independent matter haloes 
cannot produce any excess mass. On the other hand, they still may generate a GGL signal (12) if 0. If lenses cluster, i.e., a> 
or Q 0, we will get automatically non-vanishing contributions to 0. For lens-lens separations comparable or smaller than the 
typical angular size of a cluster, most lens pairs will inhabit clusters and contribute mostly to the excess mass. Therefore, those pairs 
probe essentially the matter environment of clusters, provided they are at similar redshift. We do not expect relevant contribution to 
§ from pairs of lenses with distinct redshifts (apparent pairs), though. Imagine a catalogue of lenses in which all galaxies are clearly 
separated in radial distance. On the sky, these lenses are (a) unclustered and (b) their matter environments are mutually statistically 
independent owing to the large physical distances between lenses. This exactly covers the aforementioned situation as reflected in 
a ILM with vanishing o) and Q: from this lens catalogue vanishes. In comparison with a sample of galaxies all at similar radial 
distance, a survey with radial spread in the lens distribution will have a larger fraction of apparent pairs, a reduced angular clustering 
of lenses and hence a overall suppressed amplitude of 0. This suppression can be corrected for, if the radial distribution of lenses is 
specified (Appendix A). 

The ILM is only an approximation for the clustering of galaxies and matter since a change of the matter-to-light ratio with size 
of structures is not possible. Every structure can only be a sum of individual haloes with no correlation to each other. Contrary 
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Fig. 3. Right panels: Excess mass around lens pairs with fixed separation; squares indicate the lens positions inside the maps. Left panels: Excerpts 
of the underlying ILM mock data: Lenses, shown as little dots, are either clustered (top left) or randomly distributed on the sky (bottom left). For 
simplicity, every lens has the same individual matter halo with a Gaussian lensing convergence profile (r.m.s. size is 20 pixel) sticked to it. The 
clustered lens haloes produces the joint matter halo of galaxy clusters in this model. The intensity scale in the left panels depicts the combined 
lensing convergence of all lenses; this is probed as shear by a sample of source galaxies. Note that the angular scale or the shear amplitude are of 
no particular interest here. The bottom right panel is the actual measurement of the bottom left scenario with the colour scale of the top right panel. 



to GGL this limitation is relevant for G3L, as discussed in the following. We first notice that in comparison with GGL, the lens- 
lens-shear correlator under ILM assumptions seemingly does not provide any fundamentally new information about the lens-matter 
connection. To clarify this, the clustering coiTelation functions (w, Q) can be determined by the observable lens angular distribution 
without exploiting the gravitational lensing effect. Then, by utilising u and the observed mean tangential shear y,, the average halo 
shear profile y|, can be constrained from Eq.(12) so that all essential ingredients for predicting Q, Eq. (20), are already fixed. In 
particular, Q appears to be only sensitive to the average halo shear profile as GGL is. Seen this way, G3L can at most complement 
constraints on the mean shear profile y^. In a scenario more complex than the ILM, however, this differs. Imagine throwing in a 
few completely different matter haloes into a ILM, hosting several galaxies simultaneously, that cannot be described as sums of 
individual haloes. Conventional GGL would be unable to detect a difference to a ILM scenario, as we can still use the clustering 
correlation function w and the GGL signal ji with (12) to define an average lens shear profile %. Therefore the ILM is always, 
even when falsely presumed, sufficient to consistently describe the GGL signal and the lens clustering. On the other hand, when 
then combined with Q, we would observe inconsistencies, as we fail to coiTectly explain Q with Eq. (20). In fact, the new haloes 
with more than one galaxy would produce a one-halo term (Appendix B), which is missing in the ILM description. From this we 
conclude that Q, in combination with GGL, enables us to detect whether a ILM sufficiently explains the data or whether a more 
advanced description is required. 

With the ILM as reference scenario we suggest to construct a test for the applicability of the ILM with the excess signal AQ 
constructed as follows: 

1. Use GGL, lens clustering statistics and numbers to obtain the ILM parameter set («g, w, Q., y^) via Eq. (12); 

2. Define the ILM excess signal by 

AQi&l , ^1, 03) := 0{-»i , 1?2, <^3) - Qi'^{-»x,^2, h) - Sihi&l , '?2, 03) , (26) 
where the last two terms on the rh.s. are the two- and three-halo term, Eqs. (21) and (23), from the ILM description. 
A vanishing tests the validity of a ILM description for the data or expresses the deviation from it. 
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4. Shear-shear-lens correlator 

Here we predict a measurement for the second G3L correlator G±, given through the correlation of two shears and one lens number 
density 

(rc(ei)rc (»2)«g(03)) = h-^'e^^'^'^'^'^'-^Gd^U^l, 03) . (27) 

The geometry of the correlator is depicted in the middle panel of Fig. 2. Here and in the following equations, a superscript "±" 
as in means -y^ for and -y^ for the complex conjugate y*. This correlator measures the shear-shear correlations as function 
of lens separation. As before with 0, symmetries demand that the correlator depends only on relative separations and angles given 
by the triangle defined by lens and source positions. It contains an unconnected part that describes the shear-shear correlations for 
randomly distributed lenses with no correlation to the shear field, namely (Schneider & Watts 2005): 

Gl'{i}u-»2,cf>3) ^+()?3)e+-'^-^ ; G"^(^i,i?2,03) ^-(^^3) (^e^'^'/^ _ ^e-'^'/^j . (28) 

j?3 denotes the source-source separation. As before, the unconnected terms do not contribute to the aperture statistics, here (TVM^p), 
and are thus of no particular interest for G3L. Subtracting the unconnected terms leaves us with the relevant excess shear-shear 
correlations about lenses, formally (cf. Eq. 1) 

0^(^1,^2,03) = (r(ei; ^1)7^(^2; ^2)^8(^3)) = Gi(l?i,1?2,03) -G^^()?i,1?2,03) . (29) 



4.1. Derivation 

The evaluation of G± for the ILM boils down to evaluating the connected terms of the triple correlator 



(rc(ei)re(^2)«g(03)| (30) 

= Z (^D^(^3 - e^)7m - e^-; a,)r^(02 - (^■, ai)) + 2 (5^^\03 - <)rh(0i - e)- aj)y^ie2 - e\; a,)\ 

' ^ i^J^k 

one-halo three-halo 

+ J] {s'^\03 - e))y^(ei - e); aj)y^(e2 - 0); ff,)) + ('5d'(»3 - 0))yb(0i - e); adyt{02 

two-halo-2 



+ Vo(^^ - ^5)^(61 - e^ adyti02 - e); a,)) , 



two-halo— 1 

which now contains a one-halo term, two-halo terms and a three-halo term. We distinguish two categories of two-halo terms: In 
"two-halo-I", the two shear signals are associated with the same halo, while in "two-halo-2" the shear signals originate from the 
lens halo and a different neighbouring halo. In analogy to the calculations for the correlator Q, only terms associated with the leading 
order clustering correlation functions in the halo correlators (. . (terms with cS) and (. . (terms with Q) are of interest for 
the connected terms; for (. . .), all terms are connected. Going through the averages step by step and collecting the connected terms, 
yields as final result for 0\ + Oni^ O3: 

GA&U&2, 03) ^Gi\&u&2,4>3) + Gi" (&i ,&2,4>3) + Gf\§u &2, 03) + Gf{l»i , t?2, 03) ■ (31) 

We start with the three-halo term, which is after performing the integral variable transformations 9 t-^ + 63 and 0' 6' + 63: 

G'^{&u&2,<f>3) := H2e-2'(>^'^'^=> J d^0d^0'n{\0 - 03\,\0' - 03\,\0 - 0'\)n(.Oi - O)7ti02 - 0') (32) 

^•^^^ J d~0d^6'n{\0\, m, \0' - 0\)%(0i3 + ^)rh(»23 + 0') . 



-2^-2i{ip,Tif 



To cast this into a form that no longer explicitly contains any yp,, we need to do a few more transformations. We first note that for a 
shifted tangential shear one has: 

+ = - (^7^ ^'."(1^ + 0'\) = -e"-'^'A(^, ^ - rt,h('r(0, 6', ip - ip')) , (33) 
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where 

ie^"^ + 1 1 + 2 JCOS0 + s^cos(20) - i2i(l + icosrfi) sin0 .» , . i , 

A(5,0) - —r-T = 7— ): ; A*(5,0) = A-'(5,0) (34) 

ie + 1 1 + 5-^+25 cos 4> 

is an additional phase factor; and (/?' are the polar angles of and 0' , respectively; T is given by the previous Eq. (24). Then this 
allows us to rewrite the previous equation for G^^ as 

GiVl,'?2,03) 

ji^^^Hm^ j d^0d^0'n{\0\A0'\A0' -0\)^A{j-,<f-Vi^y,^MOi3 + O\)$^ 

d99 de'e'J^ d^ d^'Al— ,^jA^(— ,^Jn(0,0\m0\^'-^ + 03))r,,(T(7?i,0,^))r,,h(m (35) 

The (connected) two-halo terms can be split into two sub-groups. One group is insensitive to shape variations of the lens halo, 
as it only contains the mean halo shear profile y^, 

Gf\^i,'9i,4>2) Hge-2i('^'^'^^' J d^9 (^[^ - 0\){^^{ei - ei)yt{02 - 0) + yk{ei - 0)yi{02 - 0^)) 
- Hge-''(^''^^' J d^9a>{e)i^^{0,,)yi(02i + 0) + r^(^23)rh(»i3 + 0)) 

Xoo pin 
d09y,,^(9) d^cos(2^)w(T(^2,0,^)) 

deOyiMiS) d^cos(2^)w(Y(^i,0,^)) 

= r,,h('?i)(rt('?2) - r,.h('?2)) + r,.h('?2)(rt('?i) - AhC'^O) ■ (36) 

The last step exploits the two-halo term of our previous result (12) for GGL. The terms inside the brackets express the excess 
tangential shear due to haloes clustering around the lens halo. The steps in the calculation of G|''~ are by and large identical to the 
steps undertaken in Sect. 2.2. 

The remaining terms in (31) are special, as they are indeed sensitive to halo variations, which sets them clearly apart from all 
aforementioned correlation functions. The (connected) one-halo term is 



Gi\&u&2, h) e-2'(^.^^2)^^^(0, _ _ Q^. ^)^^ ^ e-2'(^' ^^^'^7,(013; a)y^{023; ff))^ . (37) 

If we write the halo shear as sum of the mean shear profile and some fluctuation 6yh{0', a) about it, i.e., 

rh(0;a) = rhW + <5rh(e;a), (38) 

with {5yh{0', a)}a - 0, then the one-halo term becomes 

Gi\)?i, i?2, <p3) = e-^'t^-^-^^) Jy^(0j3)^(023) + I^Sym^; a)6y^i023; a)^j 

= rt,h('?i)rt,h('?2) + e-2'('^'^'^=)(5rh(ei3; a)Syi(023; a)) 



= r,h('?i)rt,h('?2) + {Sy,,ui&u<fiua)6yf^^{&2, n; a))^ ■ (39) 
In the last equation, we employed 

yt,h(0, ip;a):^ -c'^f y^{0- a) ; y, ^(9) = U^{9, ^; a)) , (40) 



where is the polar angle of 0\ an equivalent definition is employed for the fluctuations (5yt,h- Owing to statistical isotropy of the 
shear field, the model halo shear profile yb{0\ a) has to have a random orientation. Therefore, the correlator in the previous equation 
must be invariant with respect to any rotation 0, or 



" '^'^(^^'.h('^i'^ + '^3;a)5rth('^2,<^;a))^ 



rt,h('?i)rt,h('?2) 5Gi"(i?i, i?2, 03) . (41) 
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The remaining two-halo term in (31) is similar to the one-halo term, actually an integral over G^_^: 



J d^0co(e) ^yh(0i3 + 0; a)y^(023 + 0; a)j^ 



where the angle v spanned by 023 + and 613 + and the corresponding angle fi for a rotated by ^2 are implicitly given by 

_ 023 + \0n+0\ _ »?2e-'^^-^^> + g T()?i.g,y-yi) _ _ t^ze-'-^ + e Y()?i,e,y + ^3) 
^ 6»i3 16*23 + 6*1 )?ie-i('^-'^i) H-fii T(i?2,6i,^-(^2) ' ^ §ie-^(f+^^^ + 6 T(i?2,6',^) 



4.2. Interpretation 

The resulting G± is the lowest-order galaxy-galaxy lensing correlation function that, at least within the framework of the ILM, 
is sensitive to variations among shear profiles of haloes. Nevertheless, we also have a G± signal when all halo shear profiles are 
identical, i.e., ^yt.h = 0, which is generated by the tangential shear of the lens halo and the excess tangential shear of clustering 
neighbouring haloes. Only (39) and (42), which both contain the correlator 60^^, are affected by a scatter in halo matter profiles. The 
one-halo term of G± remains unchanged, even if we allow for general correlations between halo parameters a of distinct haloes or 
for correlations between lens positions and a, since a statistical independence of haloes has not been used for this term. Therefore, 
we expect the behaviour of G± on small angular scales to be described generally by G^^, not just within the ILM. 
In the simplest case that all halo shear profiles are exactly identical, yh{0', a) = y\^{0), we find 

Gi\§i,&2, h) = rt,h('?i)rt,h('?2) ; sg'\&u&2, (44) 

i.e., the one-halo term has no explicit dependence on the opening angle ^3. Note that G^_^ and Gi*" are identical. As illustration of 
the impact of variance in halo shear profiles, consider a singular isothermal ellipsoid (SIE) profile with ellipticity eh and random 
orientation (p (Mandelbaum et al. 2006a) 

rr(«) (1 + y + 20)) , (45) 

(p is the polar angle of 0. The absolute amplitude of the shear profile is not of interest here. The tangential shear profile of the SIE is 

rtt(^' ^) ^ (1 + I + 2^^)) . (46) 

In this case, we find by marginalising over all orientations <p: 

Gi\)?i,l?2,<^3) = ^ j| d0rt"(,?l,03)[r,t('?2,O)f (xy, ^^^^ ^'^ .- . , ^ - , ^, 1 

thus similar to the previous result but now with an additional 03-dependent term, or 



l + ^cos(203)| ;rt,hW<^-, (47) 



,5GlVi,)?2,03) = v'^os(203)- (48) 
o 

The result becomes somewhat more complicated for general slopes 5 (Appendix C) and will reveal a difference between G^^ and 
Gi^ when 5 I (not SIE) and eh it (elliptical). Moreover, the correlator (41) will exhibit no (^3 -dependence, if the lens haloes 
are always axially symmetric, even though their radial matter density profile or their mass may scatter as to be expected in reality. 
Therefore, we conclude that G± may in principle be used to constrain the shape or, more specifically, the mean second-moment of 
the projected halo matter density profiles. In addition to that, fluctuations ^yt.h in the halo shear profile due to halo substructure also 
add to the variance dependent one- and two-halo term of G+. As with the foregoing Q, it may be useful to define an excess AG±, 
obtained by subtracting off the G±-signal as anticipated from the ILM with parameters from lens clustering and GGL. In the ILM 
regime, the excess signal AG± exactly vanishes, if there is no scatter in the (projected) matter density profiles. 

5. Conclusions 

In order to gain a better understanding of G3L, we conceived a toy model, the "isolated lens model" (ILM), for the distribution of 
galaxies and matter about galaxies. In this picture, "isolated" galaxies are surrounded by their own matter envelope (halo). Variations 
in the halo matter density profile are explicitly allowed, albeit statistically independent to variations of other matter envelopes or 
to positions of other lenses. Consequently, the matter environment of clusters is herein the superposition of independent haloes 
produced by clustering galaxies. The average independent matter halo is described by the mean tangential shear around lenses and 
the clustering of the lenses (GGL), Eq. (12). The foregoing calculations evaluate what would be measured by G3L (Eqs. 20 and 31) 
under the ILM assumptions and discuss the results. Here we summarise our main conclusions. 
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5. 1 . Excess shear about lens pairs 

The lens-lens-shear correlation function basically stacks the shear field about clustered lens pairs as opposed to GGL, which 
stacks the shear field about individual lenses. "Clustered lens pairs" refers to the fact that the connected part of the lens-lens-shear 
correlator does not include the expected shear pattern around pairs formed by randomly distributed lenses.* "Excess mass" refers 
to the convergence map that corresponds to the excess shear of lens pairs. Our conclusions are: 

- Unclustered lenses do not generate any signal, although they may exhibit a GGL signal. ^ is a probe for the matter environment 
of clusters (or groups), probed by lens pairs inhabiting the cluster. 

- Apparent lens pairs formed by lenses well separated in redshift, overall diminish the signal and add noise. The signal-to-noise 
of can thus probably be improved by exploiting lens redshift information and by giving more weight to lenses that are close 
in redshift. The susceptibility of the amplitude as to the radial distribution of lenses can be calibrated (Appendix A). 

- In the ILM, the excess mass constituents are the haloes about the clustered lens pairs (two-halo term; Eq. 21) and a third halo 
clustering about the clustered lens pair (three-halo term; Eq. 23). In a more elaborate model incorporating haloes hosting more 
than one galaxy at a time, we expect also one-halo terms adding to the excess mass (Appendix B), especially at small scales. 

- If the ILM is a fair description, then does not provide any new information on the galaxy-matter connection compared to GGL 
combined with 2"''-order galaxy clustering. can at most complement the information on the mean halo shear profile. 

- G3L in combination with GGL and lens clustering can probe whether the matter environments of clusters could be expanded as 
sums of independent haloes, as in the ILM. In particular, a sum of independent haloes would be unable to change the mass-to- 
light ratio compared to that of field galaxies. 

- Moreover, the ILM can be employed as reference to quantify the deviation from the ILM picture in the real matter distribution 
around galaxies. This can be devised as practical test for a more advanced halo-model picture that naturally presumes the 
possibility of genuine joint matter haloes of galaxies that are fundamentally different to sums of matter haloes about isolated 
lenses. In analogy to 0, ISQ can be visualised as mass map about lens pairs or as aperture statistics. 

5.2. Excess shear-shear correlations about lenses 

UnUke GGL and 0, which both stack shear about lenses, G± measures the excess shear-shear correlations relative to a lens position. 
"Excess" means in this context shear-shear correlations in contrast to randomly distributed lenses, which would have a vanishing 
GGL or signal. In this respect, G± more resembles the shear-shear correlation function ^± utilised in cosmic shear studies (e.g. 
Schneider 2006) rather than GGL. Our conclusions are: 

- In the ILM, or more generally in a full halo model picture, the excess shear-shear correlations have two basic contributors: (i) 
the host halo of a lens (one-halo term; Eq. 39), and (ii) neighbouring haloes clustering about the lens (two-halo terms: Eqs. 42 
and 36; three-halo term: Eq. 35). 

- G± is the lowest-order galaxy-galaxy lensing correlation function that is sensitive to variations in the (projected) density profiles 
of matter around lenses. Traditional GGL and are only functions of the mean, stacked shear profiles. 

- In particular is G± sensitive to variations due to elliptical haloes with random orientations. Therefore, the correlator is principally 
sensitive to the shape of matter haloes and could be exploited as such to measure halo shapes without the need of luminous tracers 
of presumed alignment to the halo. An elliptical halo generates an extra 03-modulation in G+, where (^3 is the opening angle 
between the two lens-source directions. The one-halo term of G± does not exhibit a 03-modulation for spherical haloes. 

- Halo substructure, i.e., random fluctuations about a smooth halo profile, also contribute to the one- and two-halo term of G±. 
Therefore, G± is in principle also sensitive to halo substructure. 

- The one-halo term of G± is unchanged, if we generally allow for statistical dependences of the lens haloes within the framework 
of a general model. Therefore, the ellipticity and substructure effect will also be present to some extend, if we have strong 
deviations from the ILM assumptions. However, it is unclear at this point how strong the effects are, even within the ILM, and 
what possible degeneracies are. We defer a thorough study of these effects to a future paper. 
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Appendix A: Calibration scheme 

Let 

be the critical surface matter density for lenses at comoving distance ;ifd and sources at^s, and the critical density of the Universe, 
respectively; /kO^) denotes the angular diameter distance. Using Limber's equation for projecting the 3D bispectrum to the angular 
2D bispectrum, Schneider & Watts (2005) showed for 

\k\^ := ki+kl-2kik2 cos (A ; k\R\+klRl-2kik2RiR2 cos (1/^3 - i/r) ; e"'" := ^ [2kik2RiR2 + + (;t2-R2)V*"'*'](A.2) 

that 

n n r AxAx'p^{x')Pi(xf r dk,k,dk2k2d^' jki^-'^'' + k2e^'^"f ^ ^ ^ , \ ,a 
Jo I-crnix ,X)a(x) J (2^) , \kr , ^ ' 

= :K{ki.ki,Rl,R2-ili) 

Pf(x) and ph(x) are the radial lens and source distribution; Bggin{ki,k2,ifr;x) is the galaxy-galaxy-matter bispectrum at comoving 
radial distance x\ the wave numbers A:, are also in comoving units. Note that we here, inside the integral, transformed angular 
separations to the projected comoving distance at lens plane distance^, Rj :- fK(x)'^i- We recast this equation into 

a A.\ o r" <ixdx'pb(x')Pf(xf ^, fd^iM^aM^^,, , p p ,^ \ 
g{RuR2,(p3) := i^mPcrit — — —; — , , , X K{kuk2,RuR2,tl^)BggJkuk2,t//-Xi) (A.4) 

/dkikidk2k2ddf , \ 
—^^m^ ,k2,RuR2,^) BggJki , ^2, <A;^d) 

where Bggm(^i, ^2, ^'^Xdj the bispectrum at effective lens plane distance defined by 
p dxdx'pi,<x')Pf(j(f , . , \ R t , - \ r dxdx'Ph(j^')PfM' 

-^^-rr-rT^^sg4ki^k2,tf^-,x)-- Bggjki,k2,ii'-,Xd) \ ^ , , , , ^ ^^-^^ 

and R, - f^Sjld'^i is the projected comoving distance of angular separation at the effective lens plane distance. Therefore, 
normalising Q by ^0 corrects the correlator for the amplitude reduction due to lens pairs with lenses at distinct redshifts, encoded 
in pfCf), and removes lensing related quantities, yielding a bispectrum at an effective distance projected onto the lens plane by 
kernel /:(...). 
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Similarly, we find for the shear-shear-lens correlator: 

G^{Ri,R2,<p3)^G^,x ^^t}h}^^K^i^kiRuk2R2,^)B^^J^kuk2,<p\XA) (A.6) 
with 

Jo 2;critO',^)2flCr)- Jo ^cvitix ,xra(x)- / °v / ^iplnt 

and the integral kernels: 

K+[kiRuk2R2, iA) e^^'^'^-'^'^^VoCA) ; K^[kiRuk2R2, if') := e^'^'^JM) . (A.8) 

Appendix B: One-halo term for ^ in a more complex model 

The ILM does not have any one-halo terms for since haloes are hosting only one lens at a time. Matter haloes of clusters have to 
be represented as sums of independent haloes. In a more realistic situation beyond the ILM, we can imagine genuinely new haloes, 
hosting more than one galaxy, that cannot be expanded as sum of independent haloes. For simplicity, for this new type of haloes we 
here assume round haloes with shear profile jf^iO); 6 - Ois the centre of the halo. Every halo hosts the same number of A^g galaxies. 
Galaxies belonging to the /th halo are scattered throughout the halo with a separation \6ij relative to the halo centre 6^. Now the 
number density of galaxies on the sky and the (relevant) shear field associated with the new halo type are 

"g(^) = Z 2 " ^' " ^^'■'■^ ' ^^^^^ = Z - ^) ■ (B.l) 

/=1 j=l /=1 

In the following, we will focus on the connected one-halo terms of the lens-lens-shear correlator only: 

{n^{0i)n^(02)yc{0i)) = Z ("^D^^i " " ^0ij)6u{e2 - O) - ^eij)j^{ei - ^,)). . (B.2) 

'■=1 J=l '"' 

Nh 

+ YjYj (•^D^^i " ~ mj)6Y,{02 Afe)rh(03 - ei)) (B.3) 

+ C2-halo terms) + C3-halo terms) , (B.4) 
where we have the ensemble averages 

^ . . . y := ^ J d^Qi d-A0,v P^{mj\)[ ■ ■ ■] ; ( . . . ) ^ := ^ J d^fl,. d^A^?,^ d^Aft* P^{mjV mk\, mj - Mi,\)[ . . .] . (B.5) 

By Pa(A0) we denote the p.d.f. of a single relative galaxy position inside a halo, whereas by P/^{A6, AO', \A.6 - A.6'\) we denote the 
joint p.d.f. of two galaxy positions. Owing to isotropy, the former can only be a function of the modulus of AO, the latter only a 
function of the relative separations |A^|, \A6'\ and \A0 - Ad'\. The mean number density of galaxies populating the new haloes is 
NhNg/A. If the new halo type is mixed with other haloes, such as ILM haloes, the total number density of galaxies, «„, may be 
different. The one halo term is not affected by the presence of the other haloes apart from /Ig. 

The sum (B.2) vanishes for 9i + 62 so that the only relevant contribution to the one-halo terms of is, after exploiting the Dirac 
delta functions: 

@'\M2,4>3) = -e-'<'^'+'^^'«gXA^g(A^g - i)A-' J d-ep^[\e + en\,\e\,\0n^MOn + e) 

= -e-'('^'+'^%2g7hA^g(A^g - 1) J d^0PA{\0 + 0nl\0 + 023l&^)%{0) 
= «g7hA^g(A^g - 1) J d0W^PA(|0 + ^i3|,|^ + 023l,'?3)e"2^"^>-^^'rt,h(e) 

poo p27r 

'fhNg(Ng-l)J^ d06%^,(0) dip cos(2<p)PA[r(&uO,<p + (f>i),r(§2,0,<p),^i) . (B.6) 



We used (24) for the definition of T. By /h :- Nh/iAiig) we mean the ratio of new type haloes to the total number of galaxies; 
as usual, = -1- )?2 - 2d-\§2 cos ^3 is the separation of the lenses. If we add more haloes with, say, different shear profiles or 
numbers of galaxies A^g, we will obtain a sum of one-halo terms of the previous kind, all weighed with (i) their halo fractions /h 
and (ii) number of galaxy pairs Ng{Ng - 1). ILM haloes have trivially A^g = 1, thus vanishing one-halo terms due to the absence 
of galaxy pairs. Moreover, if galaxies are unclustered inside their host haloes, P^ ~ const, one will also have a vanishing one-halo 
term. When the host haloes are not clustered as well, there will be no contribution to at all. 
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Appendix C: General elliptical matter haloes and G+ 

We consider a projected halo matter density profile with constant power law slope when averaged over annuli. The profiles have an 
elliptical shape with ellipticity - (a^ - b^)/{a^ + b^); a,b are the sizes of the major and minor axis, respectively. Mandelbaum 
et al. (2006a) find for the shear profile corresponding to this matter density profile 

rh(0;a) = -— e ^(y. + iy,); 1 + COS (2y) ; yx :^ Ysm(2y), (C.l) 

where we denote by the polar angle of and by A the amplitude of the shear. For simplicity, we work with the assumption that 
all haloes have the same ellipticity, power law index and amplitude. As the orientation of a halo is a-priori not known, we have to 
marginalise over the random orientation angle of the halo for correlator (41): 

—^{B,^) = ^-7 , (C.2) 

- I 

2 

,5GiVi,'?2,03) = ^/-(5) COS (203), (C.3) 
o 

2 2 

5GlVi,'?2,03) = ^A(5)cos(203) + i^^+(5)sin(203), (C.4) 

o o 

with the auxiliary functions 

{x-lfix + T) (/-4x3 +28x2- 48x + 32)(x- 2)2 ^ (x - 2)V - Sx^ + 6x - 4) 

x-'(x-4) (x-4)^x^ (x-4)'^x^ 

One has f+{6) = f-{6) = 1 and g^(6) = for a SIE (5=1). 
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